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Band structures for solid rare gases (Ne, Ar) have been calculated using the GW approximation. All electron 
and pseudopotential ab initio calculations were performed using Gaussian orbital basis sets and the dependence 
of particle-hole gaps and electron affinities on basis set and treatment of core electrons is investigated. All 
electron GW calculations have a smaller particle-hole gap than pseudopotential GW calculations by up to 0.2 
eV. Quasiparticle electron and hole excitation energies, valence band widths and electron affinities are generally 
in very good agreement with those derived from optical absorption and photoemission measurements. 
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I. INTRODUCTION 

Optical spectra and band structures of rare gas solids (RGS) 
have been studied, both experimentally and theoretically, for 
over 40 years. Their importance lies in the simplicity of 
their crystal structure, the nearly atomic character of valence 
states versus extended character of conduction states and the 
fact that they have strong many-body effects in their optical 
spectra. They are an important testing ground for electronic 
structure methods and as electronic structure methods have 
developed, they have been applied to RGS. Early electronic 
structure studies included applications of density functional 
theory (DFT)ii2i^ Hartree-Fock theory (HFT)i^ and self- 
interaction corrected DFT— . Some of these studies have in- 
cluded correlation effects in the band structure via many-body 
perturbation theory ''^i^i^i^ . 

In this paper we present results of ab initio all elec- 
tron and pseudopotential many-body calculations of the band 
structures of solid Ne and Ar using Gaussian orbital ba- 
sis sets. Band structures are calculated using the GW 
approximation'" " and the dependence of particle-hole gaps 
and electron affinities on basis set and treatment of core elec- 
trons is investigated. In another paper— we will present re- 
sults of calculations of optical spectra of these solids using a 
Bethe-Salpeter formalism^. Calculations were performed in 
a Gaussian orbital basis using the EXCITON codeJ-A, which 
is interfaced to the CRYSTAL code'"*. Single particle wave 
functions, energy eigenvalues and matrix elements of the ex- 
change correlation potential from CRYSTAL are used by EX- 
CITON to perform GW and exciton calculations. The prin- 
cipal parameters of the band structures of RGS which have 
been obtained experimentally are the particle-hole band gap, 
Eq, the valence band width, Wv and the electron affinity, Ea- 
The particle-hole gap is the energy difference for particle and 
hole excitations at the conduction band minimum and valence 
band maximum. It has been obtained experimentally from ab- 
sorption spectrumi^ and photoemission measurement a ' ' ^ . 

The GW approximation is a many-body perturbation the- 
ory and therefore contains corrections to a simpler, single- 
particle (SP) Hamiltonian. It was originally applied to semi- 
conductors using a DFT Kohn-Sham Hamiltonian, a plane 
wave basis set and pseudopotential (PP) approximation for 
core electronsii and was found to give excellent agreement 



with experiment for relatively narrow band gap materials such 
as Si, where DFT results in an indirect band gap which under- 
estimates the experimental value by 0.7 eV. 

All electron GW calculations have been performed re- 
cently for a variety of crystalline solids, including Sii^ii2i22i^. 
The all electron indirect band gaps calculated for Si were 
underestimatedi^i22iSii by 0.2 to 0.3 eV. There has been some 
debate whether this is due to incompleteness of the basis^^ 
or explicit inclusion of the core electrons (all electron rather 
than PP approximation)i^ii2i22iSLS. In the present work we 
compare results for the RGS using both all electron and PP 
approximations for the core electrons. 

The remainder of this paper is organized as follows; in sec- 
tion II the GW formalism used here is outlined, in section III 
results of GW bandstructure calculations are compared to ex- 
periment and earlier GW calculations on Ne and Ar. Finally, 
conclusions are given in section IV. 



II. GW APPROXIMATION 

A. Quasiparticle energies 

Conceptually, single electron and hole excitations are de- 
scribed 

(1) 

The self-energy operator, 'E{r,r',E), is non-Hermitian and 
so eigenvalues, Cm, have real and imaginary parts, the real 
part being the quasiparticle energy, E^^^ and the imagi- 
nary part being related to the quasiparticle lifetime. In this 
work the self-energy operator was computed within the GW 
approximation-^ in which the self-energy operator is ob- 
tained from convolution of the non-interacting single-particle 
Green's function, Go, and the screened Coulomb interaction, 
W, 

Y.{r,r',E) = ^ J c"'"^"^ Goir,r' , E-uj)W{r,r' ,uj) du;. 

(2) 
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Go is constructed from DFT single-particle orbitals, V-'n^^ 
eigenvalues, Ef^^ , of the Kohn-Sham operator, 



' I,) — 
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(3) 



For light elements it has been found that quasiparticle am- 
plitudes are well approximated by DFT wave functionsii. 
Thus, quasiparticle energies are simply given by 



E' 
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In this case only diagonal elements of the self-energy ma- 
trix and exchange-correlation potential, T4c[fit>], are required. 
Vic['^u] is the exchange-correlation potential of the valence 
electron density, from the initial DFT calculation. Equa- 
tion is solved using a scheme given by Hybertsen and 
Louiaii. 

E?nk « Ei^k + Z„Mii\mE'SM'S)- (5) 

Zmk is the quasiparticle renormalization factor and is defined 
by 



dmtl\^c{EMil) \ 

dE \p^=K 
and the operator AE is given by 

AI](r, r', E) = S](r, r', E) ~ V.,c[n,]. 



(6) 



(7) 



where is the crystal volume and e^g/ is the inverted, sym- 
metrized RPA dielectric matrix. Numerical evaluation of W 
and E (Eqs. (|3 and ^) requires calculation and inversion 
of the dielectric matrix at many values of uj. Although such 
schemes have been carried outa^ it is both time consum- 
ing and unnecessary in the present work. Instead, we adopt 
a plasmon-pole model based on the work of von der Lin- 
den and Horsch?^, which uses the concept of dielectric band 
structure^^ to approximate the frequency dependence of the 
dielectric matrix. The model assumes that all frequency de- 
pendence is projected onto eigenvalues of the inverted dielec- 
tric matrix, s^i{uj) through the approximation. 
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(11) 

Zqi are pole strengths, uigi are plasmon frequencies and O"*" 
is a positive infinitesimal. Eigenvalues of the inverted dielec- 
tric matrix determine the pole strengths and a plot of their 
dispersion with wave vector is known as the dielectric band 
structure^^; the dielectric band structure for fee Ar was re- 
ported previously- - . Plasmon pole frequencies are calculated 
either using the Johnson sum rula*^ or by fitting the dielectric 
matrix at zero frequency and a finite, imaginary frequency, 
iLUf. It can easily be shown that, 



^ql 



(12) 



B. Self-energy matrix elements 



and that. 



The screened Coulomb interaction is computed from the 
dielectric function and the bare Coulomb interaction 

W{r,r',uj)=^ J e~\r,r",uj)v{r",r')dr" (8) 

Two-point functions in a crystal lattice such as the screened 
interaction have the property, /(r + R,r' + R) = /(r, r'), 
where i? is a Bravais lattice vector, owing to translational 
symmetry. They may be represented as a Fourier transform 
asii, 

/(r,r')- e'^''+°''VGG'(9)e-*(''+°''-'-', (9) 

q,G,G' 

where G is a reciprocal lattice vector and q is a wavevector 
in the first Brillouin zone. Fourier coefficients of the screened 
potential, Wgg' are given by. 



WGG'{q,uj) 



1 



n \q + G\\q + G^ 



-e^\,,iq,u;), (10) 
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Results presented in Section Hill were obtained using plas- 
mon frequencies from the Johnson sum rule; quasiparticle 
band gaps obtained in this way were generally higher than 
those obtained by fitting by a few hundredths of an electron 
Volt. The plasmon-pole form of the inverted dielectric ma- 
trix allows the frequency integration in the calculation of the 
self-energy matrix element in Eq. (|3 to be done analytically. 
This leads to two contributions to the self-energy: an energy 
independent, Hartree-Fock exchange term. 



{mk\Yix\mk) = 

4^2e2 \lmk\e-<i+^y'-\nk + q)\^ 



q,G n 



\q + G\ 



where sum over bands, n, extends only over occupied states. 
The second, dynamic part 
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{mk\T.c{E)\mk) = 2^ zl i "t^/v t^^i i / 



q,G,G' " 

T.^L-G{-q)V*-G'i-l) 



\q + G\\q + G'\ 



2 E - Enk+q + i^-qisign{EF - Enk) 



(15) 



contains correlation energies of electron or hole 
quasiparticles^. Vi,-g are eigenvectors of the static, 
symmetrized dielectric matrix scG'iQ, '-^ = 0). 



C. Numerical details 

The starting point in our approach is to generate 
non-interacting single-particle Green's functions of an N- 
electron system. We use density functional theory^ (DFT) 
within the Perdew-Wang generalized gradient approxima- 
tion^ (PWGGA) to obtain eigenvectors and eigenvalues of 
the Kohn-Sham Hamiltonian. For this part of the calculation 
we employ the ab initio package CRYSTAljil which uses the 
Linear Combination of Atomic Orbitals (LCAO) approach to 
expand the Bloch functions. In order to investigate conver- 
gence criteria within the Gaussian orbital framework and ef- 
fects of core electrons we performed both all electron and PP 
calculations for each solid. Basis sets with 52 and 53 functions 
per atom were developed for PP and all electron calculations, 
respectively, for Ne. The PP and all electron basis sets for Ar 
contained 56 and 66 functions, respectively. Details of basis 
sets are given in AppendixlAl Pseudopotentials from Durand 
and Barthelat^ were used in PP calculations. Experimental 
lattice constants^ were used. 

The sum over q points in Eqs. J14> and (I15> as well as inte- 
gration over the Brillouin zone in the dielectric matrix calcula- 
tion is performed using Monkhorst-Pacl^ special points. The 
singularity in Eqs. \\A\ and \\5\ of l/g^ type for q ^ Q and 
G = G' = was integrated out using the auxiliary function 
technique of Gygi and Baldereschi24, while the singularity in 
Eq. M5\ of 1 /g type was neglected since the final result is not 
affected if it is neglectediS. 

Two special points in the irreducible Brillouin zone were 
used for calculation of self-energy matrix elements and an 8 x 
8x8 grid in the full Brillouin zone was used for the dielectric 
matrix calculation. Up to 400 (8000) G vectors are required 
to achieve convergence of the Hartree-Fock part of the self- 
energy (Eq. WAV ) for PP (all electron) basis sets for second 
row elements. In the summation over G and G' vectors in 
Eq. M51 . 65 vectors gave well converged results for all solids. 



D. Core-valence exchange-correlation decoupling 

When matrix elements of the AS operator in Eq. are 
evaluated, contributions from core electrons to the valence 



electron self-energy must be consideredii*ii. We compare re- 
sults from two alternative approximations for the energy inde- 
pendent part of the AE operator which were applied recently 
in all electron GW calculations on Si^^. The first approxi- 
mation is to compute matrix elements of the DFT exchange- 
correlation potential using the valence electron density only 
and to restrict the sum on occupied states in Eq. MA\ to va- 
lence states only, 

(mfclASlmfc) = (TOfc|S™'|mfc> + (TOfc|S™'|mfe> 

-{mk\Vxc[nv\\mk). (16) 

The second approximation is to replace matrix elements of the 
valence-density-only exchange correlation potential, V^ci?^?;], 
in Eqs. Q and ( I16> by 

{mk\V^c[nc + n^]\mk) - (TOfc|E™'^*=|mfc). (17) 

The notation {mk\Yj''^?^''^\mk) and (mfc|S]™'|TOfc) indicates 
that the sum on n in Eq. is limited to core or valence 
states only. Matrix elements of the LDA exchange potential 
and Hartree-Fock exchange operator for valence band maxi- 
mum and conduction band minimum states in silicon obtained 
by Arnaud and Alouani'- using a projector augmented wave 
(PAW) method and in this work using CRYSTAL are com- 
pared in Table U Remarkably good agreement was found be- 
tween LDA exchange potential matrix elements from either 
method (differences in matrix elements are only 0.01 eV in 
three out of four cases) and good agreement between Hartree- 
Fock exchange matrix elements is also obtained (within 0.1 
eV). The shortcoming of the latter approach (Eq. MIV ) is slow 
convergence of the Hartree-Fock exchange operator for core 
states, (mfc|S™''^|mfc). 

III. RESULTS AND DISCUSSION 

One of the aims of this work is to compare results of GW 
calculations on simple atomic solids which treat core electrons 
either by a PP or by explicitly including them in an all elec- 
tron calculation. Energies of states at valence band maxima 
and conduction band minima are given in Table HII as well as 
fundamental band gaps and valence band widths within DFT 
and GW approximations and experiment. GW all electron 
quasiparticle energies obtained by the two core-valence elec- 
tron decoupling methods outlined in Eqs. M6\ and Mil are 
given in columns labelled E'^-^ and E'^^ . DFT-PP calcula- 
tions underestimate experimental band gaps by 45% for Ne 
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TABLE I: Matrix elements of the Hartree-Fock exchange (Ex ) and 
LDA exchange potential (Wx) operators for self-consistent DFT wave 
functions at valence band maxima and conduction band minima for 
Si. The symbols ric+riy , core and denote whether core + valence, 
core only or valence only states are included in the operators. 

l4K + n„] SS""" Vxl^ 



Tisc" -11.75 -1.32 -10.18 

risc* -11.74 -1.40 -10.19 

Fas'/ -13.55 -1.80 -11.45 

Fas'/ -13.45 -1.85 -11.46 



"This work. 
''Referenced 



TABLE IL DFT eigenvalues and GW quasiparticle energies at va- 
lence band maxima and conduction band minima, valence band 
widths and energy gaps for Ne and Ar RGS. Calculations were per- 
formed using pseudopotentials (second and third columns) and all 
electron basis sets (fourth to sixth columns). The fifth column, E^i^ 
gives all electron GW data when l/rc[?^u] is calculated explicitly 
(Eq. il6i ) and the sixth column gives quasiparticle energies when 
Vjcf'T-u] is calculated using Eq. jl7t . The last column gives exper- 
imental data. Experimental dataTTtaken from ReferencelB 

unless 

cited differently. Energies are given in eV. 



pp 

DFT 




DFT 


all-electron 

^(1) 




Exp. 


Neon 












ri5„ -13.14 - 


-19.37 


-13.18 


-19.07 


-19.10 


-20.21 


Tic -1.35 


0.86 


-1.42 


0.97 


1.03 


1.3 


VF„ 0.71 


0.93 


0.79 


0.91 


0.93 


1.3° 


Eg 11.99 


20.23 


11.76 


20.04 


20.13 


21.51 


Argon 












Fis. -9.74 - 


43.15 


-10.27 


-13.02 


-13.00 


-13.75 


Tic -0.60 


0.72 


-0.76 


0.80 


0.81 


0.4 


Wv 1.35 


1.73 


1.32 


1.83 


1.85 


1.7" 


Eg 9.13 


13.89 


9.51 


13.82 


13.81 


14.15 


"Reference [TtI 



and by 35% for Ar. GW-W calculations show significantly 
improved agreement with experimental data in each solid; the 
band gap error is reduced to 6% in Ne and 2% in Ar; slight 
underestimation of band gaps in RGS is similar to that in 
semiconductors where, for example, the band gap is under- 
estimated by 4% in Si2i. The reason for good agreement be- 
tween quasiparticle energies, E^^ and E^^ , for Ne and Ar 
(this work) and SUS is that the energy independent part of the 
AE operator is E™' - KcK] for ^ and it is 
- Vxc[nc + ny] for i?^^; the difference in these two quanti- 
ties is of order 0. 1 eV and results in nearly equal quasiparticle 
energies, and E§^. 

There is good agreement between PP and all electron DFT 
calculations for Ne; the bottom conduction (Fis^,) and top va- 
lence band (Fic) energies and valence band widths (Wy) differ 
by less than 0.1 eV. Fi5„ and Tic quasiparticle energies from 



all electron calculations lie slightly above PP values. The ab- 
solute value of the Tic conduction band energy, which deter- 
mines the electron affinity, has the correct sign in GW calcu- 
lations and lies just 0.4 eV below the experimental value for 
the electron affinity, whereas DFT calculations predict larger 
electron affinities of the wrong sign. GW calculations result 
in valence band widths ^ 0.9 eV, which are smaller than the 
experimental value of L3 eViii, but are in agreement with the 
value of 0.99 eV obtained by Bacalis et ali The two methods 
used for core-valence decoupling (Tablellll columns 5 and 6) 
result in Fisu and Tic quasiparticle energies which differ by 
only - 0.05eV. 



The Fi5i, valence band maximum state in all electron DFT 
calculations on Ar is lower than in PP calculations by 0.53 eV 
while the Tic conduction band minimum state is lower by 0. 16 
eV. However GW quasiparticle energies for these states using 
either all electron or PP basis sets are in good agreement, the 
maximum difference being only 0.13 eV. The GW Tic con- 
duction band energy exceeds the experimental electron affin- 
ity by ~ 0.4 eV, whereas the DFT Tic energy again predicts 
an electron affinity with the wrong sign. GW valence band 
widths of 1.73 (PP) and 1.83 eV (all electron) agree very well 
with the experimental value of L7 eV— . The two methods 
used for core-valence decoupling also result in very similar 
quasiparticle energies for Ar. 



GW band structures along A and S symmetry lines for Ne 
and Ar are shown in Fig. ^ Self-energy coiTections to GW 
band structures in both Ne and Ar are relatively independent 
of wavevector, leading to a scissor type opening of the band 
gap on going from DFT to GW energy bands. DFT band- 
structures are not shown in Fig. ^for clarity. TablesHIIIand llVI 
give a direct comparison of all electron DFT energy eigenval- 
ues and GW quasiparticle energies at high symmetry points 
for Ne and Ar and include results from Bacalis et air and ex- 
periment. When DFT and quasiparticle energies for Ne at X 
and L points are compared (Table HIH. we find a widening of 
the valence bands by approximately 30%. Our results for va- 
lence band energies and widths are in very good agreement 
with those reported earlier by Bacalis et ali. A similar pat- 
tern of valence band widening for GW valence bands in Ne is 
found in Ar and our results are again in good agreement with 
those of Bacalis et al. The energy difference of the first and 
second conduction bands at the F point (F25C' — Tic) in our 
GW calculation is 2 eV smaller than in the PAW^ calcu- 
lation for Ne, while the value of 8.21 eV for Ar agrees well 
with the value of 8.44 eV obtained by Bacalis et al. . This 
energy difference is sensitive to completeness of Gaussian or- 
bital basis sets (see AppendixIXIi as the F25C' state has signif- 
icant amplitude in octahedral interstitial regions. Inclusion of 
interstitial functions in basis sets (Appendix 1X1 and optimi- 
sation of the most diffuse functions reduced the F25C' — Fic 
energy difference significantly, while basis sets with no inter- 
stitial functions result in a larger conduction band separation 
and fundamental gap. 
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FIG. 1: GW band structure for Ne (left panel) and Ar (right panel). 



TABLE III: Energy eigenvalues in eV at high-symmetry points for 
solid Ne. The reference energy is the valence band maximum en- 
ergy. Results in the second and third columns were obtained using an 
all electron basis set and valence-core electron decoupling was done 
using the method outlined in Eq. 1161 Results in the fourth and fifth 
columns are from all electron PAW calculations-. The last column 
presents experimental values. 





This work 


PAW 


a 






DPT 


GW 


DPT 


GW 


Exp. 




0.0 


0.0 


0.0 


0.0 


0.0 


Tic 


11.76 


20.04 


11.40 


16.56 


21.51* 




-0.61 


-0.82 


-0.67 


-0.88 




X5V 


-0.21 


-0.29 


-0.23 


-0.30 






-0.69 


-0.91 


-0.75 


-0.99 


-1.3^ 




-0.07 


-0.10 


-0.07 


-0.09 




^25c — Tic 


17.85 


18.57 


18.12 


20.51 




Xic — Tic 


6.78 


6.66 


6.82 


8.12 




Lie — Tic 


5.57 


5.91 


6.03 


7.21 





"Reference 3 
^Reference 1, 
'^Reference 1' 



TABLE IV: Energy eigenvalues in eV at high-symmetry points for 
solid Ar. The reference energy is the valence band maximum energy. 
Results in the second and third columns were obtained using an all 
electron basis set and valence-core electron decoupling was done us- 
ing the method outlined in Eq. 1161 Results in the fourth and fifth 
columns are from all electron PAW calculations-. The last column 
presents experimental values. 





This work 


PAW 


a 






DPT 


GW 


DPT 


GW 


Exp. 




0.0 


0.0 


0.0 


0.0 


0.0 


Fic 


9.51 


13.82 


8.09 


11.96 


14.15* 


Xav 


-1.19 


-1.65 


-1.28 


-1.73 




Xzv 


-0.42 


-0.50 


-0.46 


-0.63 




L2V 


-1.31 


-1.83 


-1.41 


-1.92 


-1.7^ 


/ 

Lsv 


-0.14 


-0.20 


-0.16 


-0.20 




T25C - Ti 


6.85 


8.21 


7.43 


8.44 




Xic — Tic 


2.02 


2.04 


2.63 


3.10 




Lie — Tic 


2.25 


2.56 


2.94 


3.50 





"Reference 3 
''Reference 1, 
'^Reference 1 



IV. CONCLUSIONS 

Band structures of solid Ne and Ar have been calculated us- 
ing the GW approximation. Calculations were performed us- 
ing experimental lattice constants. Gaussian orbital basis sets 
were used throughout and core electrons were treated either 
explicitly with all electron basis sets or by pseudopotentials. 
Results of all electron and pseudopotential calculations are in 



good agreement, although the fundamental band gap predicted 
by all electron calculations is smaller than that in pseudopo- 
tential calculations by up to 0.2 eV. Positions of conduction 
band minima for Ne and Ar in GW calculations are in good 
agreement with experimental electron affinities so that abso- 
lute positions of quasiparticle energy levels in Ne and Ar are 
reliably predicted in the GW approximation. Fundamental 
band gaps for Ne and Ar are in good agreement with experi- 
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TABLE V: Basis sets used in this work. Exponents of s, p and d 
Cartesian Gaussian orbitals which were centered on the nuclear site 
((0,0,0) labeled Nuc.) and at the octahedral interstitial site of the fee 
lattice ((0.5,0.5,0.5) labelled Oct.) are tabulated in atomic units. Ba- 
sis sets for atomic cores in all electron calculations were conventional 
quantum chemistry core basis sets and are not given here. 



Ne 
Nuc. sp 
Nuc. d 
Oct. sp 



32.0 16.0 4.0 



Ar - Basis Set 1 
Nuc. sp 32.0 
Nuc. d 
Oct. sp 
Oct. d 

Ar - Basis Set 2 
Nuc. sp 85.0 
Nuc. d 
Oct. sp 
Octd 



16.0 4.0 



2.0 



2.0 



34.0 14.0 



1.4 
1.05 



1.4 
1.6 



1.4 
0.8 



0.8 
0.79 



0.47 
0.8 



0.47 

0.4 

0.47 



0.39 
0.39 
0.61 



0.185 

0.2 

0.3 



0.15 
0.2 
0.15 
0.4 



0.2 
0.1 
0.31 
0.2 



mental gaps from photoemission and optical absorption data 
where shifts in the gap due to electron-hole attraction have 
been subtracted. 
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APPENDIX A: BASIS SETS 



The construction and use of an appropriate basis set con- 
stitutes a critical factor in ab initio calculations and is partic- 
ularly important within a Gaussian orbital framework. Apart 
from minimizing the total energy, which is necessary for a 
good quality basis set, one has to ensure that the basis set con- 
tains a sufficient number of basis functions to generate the vir- 
tual space. The need for a large number of conduction bands 
for a well converged self-energy has been emphasized again 
recently22i. The number of conduction bands can be increased 
by including / and g type functions into the basis set, but these 
are not yet available in the CRYSTAL code. Alternatively, ex- 
tra sets of orbitals were added at interstitial sites of the crystal. 
This improves the flexibility of the basis set through the unit 
cell and attempts to reproduce the highly nodal structure of 
free-electron conduction band states. 



Two techniques were used for constructing basis sets: 
Firstly, starting from two decay constants, 0.15 and 2.0, ge- 
ometrical expansion was used to generate more localized or- 
bitals, interstitial functions were added and the most diffuse 
functions were adjusted to minimize the total energy. The 
second approach used valence exponents from conventional, 
contracted quantum chemistry basis sets. Several Gaussian 
functions are combined into a single basis function in a con- 
tracted basis function by fixing their weights. Here the same 
exponents as used in contracted basis functions were used, but 
relative weights of different exponents were determined dur- 
ing the self-consistent field DFT calculation. The basis set 
used for PP Ne and Ar (Basis set 1) and all electron Ne calcu- 
lations was of the first type while all electron Ar calculations 
were performed using a basis set of the second type (Basis set 
2). 
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